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(D . Abstract 

We study the associated production of the neutral CP-odd Higgs boson with 
a neutral gauge boson Z in high energy e~^e~ collisions at the one loop level. We 
present a detailed discussion for the total cross-section predicted in the context 
' of the Minimal Supersymmetric Standard Model (MSSM) and make a comparison 

with the non-SUSY Two Higgs Doublet Model (THDM). We show that the MSSM 
cross-section may be enhanced compared to that for the THDM. 
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1. Introduction 



The Standard Model (SM) of electroweak interactions [|l| is in complete agreement with 
all precision experimental data (LEP, Tevatron, SLD). So far undiscovered is the SM 
Higgs boson (0°), which is the major goal of present and future colliders. Global SM fits 
[0] favour a relatively light 0°, suggesting that it should be observed soon at forthcoming 
colliders. A light Higgs boson is also a feature of Supersymmetric (SUSY) theories, which 
are the most popular extensions of the SM. Since 1989 LEP has searched for 0*^ in the 
energy range 91 GeV < ^/s < 208 GeV via the mechanism e~^e~ — *■ Z*^*)0°. In the final 
run at energies around y/s = 208 GeV two of the four experiments found excesses in the 
search for Higgs boson production in association with a Z boson (Higgsstrahlung), 
while the other two presented improved lower limits on its mass 0] of M^o > 109.7 GeV 
and > 114.3 GeV. 

LEP agreed that a further run with about 200 pb~^ per experiment at ^/s = 208.2 
GeV would be enable the combined data from the four experiments to establish a 5a 
discovery However the extended run was not approved and LEP was consequently 
shut down. The Higgs search will continue with Run H at the Tevatron which has a 
chance of confirming the existence of the Higgs boson in the mass range hinted at by LEP 
(M^o f» 115 GeV). This region is fairly problematic for the LHC and would require 
several years searching for 0° 77 decays in order to confirm such a light Higgs. 

The still hypothetical Higgs sector of the Standard Model (SM) can be enlarged and 
some simple extensions such as the Minimal Supersymmetric Standard Model (MSSM) 
and Two Higgs Doublet Model (THDM) P, Q are under much intensive study. Both the 
THDM and MSSM introduce 2 Higgs doublets which break the electroweak symmetry [^. 
From the 8 degrees of freedom initially present in the 2 Higgs doublets, 3 correspond to 
the masses of the longitudinal gauge bosons, leaving 5 degrees of freedom which manifest 
themselves as 5 physical Higgs particles (2 charged Higgs H^, 2 CP-even H^, and one 
CP-odd A^). Aside from the charged Higgs sector, the neutral scalars of the THDM may 
possess a very different phenomenology to that of the SM Higgs. The discovery of a CP- 
odd Higgs boson or/and charged Higgs boson would be clear evidence of physics beyond 
the SM. Detection of a CP-even neutral Higgs with couplings differing significantly from 



those expected for 0" would also provide evidence for physics beyond the SM. If the 



Higgsstrahlung cross-section is suppressed then the mass bound on such a Higgs boson is 
much weaker than the above bounds on M^o. 

Until now, no Higgs boson has been discovered, and from negative searches one can 
derive direct and indirect bounds on their masses. The combined null-searches from all 
four CERN LEP collaborations derive the lower limits Mh± > 77.3 GeV and Ma > 90 
GeV [|, 0] in the context of the MSSM. 

CP-odd Higgs bosons can be produced at e'^e~ colliders |14] via e+e^ h^A^ and 
e~^e^ bbA^,ttA^ [|T^,||T6[. At future e"*"e~ colliders the simplest way to produce CP- 



even Higgs scalars is in the Higgsstrahlung process e+e ^ Z* ^ HZ. The CP-odd 
A^ possesses no tree-level coupling A^T/L^ and the other tree-level diagrams for such a 
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process are proportional to the electron mass and consequently negligible [|. The dominant 
contribution is therefore from higher order diagrams which will be mediated by both SM 
and non-SM particles. Therefore the rates are expected to be strongly model dependent. 

Previous work on the loop induced production of Higgs bosons in association with 
gauge bosons (e.g. e'^e~ — > A°7, H'^W'^) can be found in 0, |1^] and references 
therein. Recently, associated production 77 — > ZH, both for SM Higgs and MSSM Higgs, 
has been studied in In a previous paper [|T^] we studied the process e~^e~ A'^Z, in 
the context of the THDM. In this paper we extend the study to the case of the MSSM. 
The particle spectrum of the MSSM doubles that of the THDM, but the scalar sectors 
are equivalent up to some mass relations and differences in the couplings 0. One would 
expect differences in the rates for the THDM and the MSSM, since in the latter there will 
be contributions from loops involving sparticles. If e"'"e~ ^ Z* ^ A^Z were sizeable it 
would provide an alternative way of producing A^ at e+e^ colliders, with a kinematical 
reach superior to that for the mechanism e~^e~ —>■ Z* ^ A^H^. Note that in the MSSM 
the kinematically favoured mechanism e~^e~ Z* ^ Al*}^ is suppressed by the factor 
cos^(/3 — a) which is very small in the region Ma > 200 GeV. The study of the various 
production mechanisms of the CP-odd Higgs boson is well motivated since the discovery 
of such a particle would signify that the electroweak symmetry breaking is introduced by 
more that one Higgs doublet (0° is CP-even). 

Our work is organized as follows. In section 2 we outline our approach for evaluating 
the 1-loop rate for e^e' A^Z in the MSSM. In section 3 we present our numerical 
results and section 4 contains our conclusions. 



2. One— Loop Corrections 

In the limit of vanishing electron mass, the process e+e~ — >■ ZA^ posseses no tree-level 
contribution and is thus mediated by higher order diagrams. We have evaluated the 
one-loop amplitude in the 'tHooft-Feynman gauge. This amplitude contains ultraviolet 



divergences (UV) and we will use the dimensional regularization scheme |^2|] to deal with 
them. 

The typical Feynman diagrams for the virtual corrections of order are drawn in Fig.l. 
These comprise: 

(i) The fermionic contribution (s-channel) to the vertices 'y-Z-A^ and Z-Z-A^: Fig. 1.1 

1.2 (+ crossed counterpart) and Fig. 1.3 

(ii) SUSY t-channel correction to the vertex e~e'*A^: Fig. 1.4, Fig.l. 5 (+ crossed coun- 
terparts). This kind of vertex is non- vanishing in the limit of vanishing electron 
mass because one of the electrons is off-shell. 

(iii) THDM box contributions. Fig. 1.6, and MSSM box contributions. Fig. 1.7 — * 1.8, 
(+ crossed counterparts) and topology Fig. 1.9 and Fig. 1.10. Note that one-loop 



^Note that at a muon collider |22l, the tree-level diagrams cannot be discarded anymore and higher 
order diagrams would induce corrections to the tree-level rate. 
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Figure 1: Generic Feynman diagrams (vertices and boxes) contributing to e~^e — > ZA^ 
in the MSSM 

contributions coming from initial state e^e'Ho, e~^e~ho vertices also vanish in the 
limit TUe 0, since e"*" and e~ are both on-shell. 

The t-channel vertex (Fig. 1.5) and box diagrams Fig. 1.8 and Fig.1.9 with exchange of 
charginos/sneutrinos contain vertices which violate fermion-number and we use Denner 
rules to deal with them. All the Feynman diagrams are generated and computed 
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using FeynArts and FeynCalc ^4] packages. We also use the fortran FF-package |^ 
the numerical analysis. 

In the MSSM, the s-channel topology in (i) also receives contributions from sfermion 
triangular loops. However, the coupling of to a pair of sfermions satisfies the following 
relation: {A^ fif*)i^j = —{A^fjf*)i-^j and consequently the set of sfermions contributions 
to 7-Z-y4° and Z-Z-A^ vanishes. 

Diagrams which involve Z-A^ (resp Z-G°) mixing in the s-channel self-energy and external 
lines vanish because the vertices j-A^-A'^ (resp '-/-G^-A^) and Z-A'^-A^ (resp Z-G^-A^) are 
absent. Z-A° mixing has to be considered in the t-channel, but owing to Lorentz invariance 
the Z-A^ self-energy is proportional to CP odd momentum p'^o = {Pe+ + Pe- ~ VzY] 
then, since the vector boson Z is on-shell, the t-channel amplitude will be proportional 
to irif. and consequently vanishes. Note also that the tadpole diagrams have vanishing 
contributions. We stress in passing that, there are diagrams in the process e^e~ A^Z 
which do not contribute to e^e~ A^'^. These are diagrams which involve the couplings 
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Zx'^X^ S'lid Zui), which would be absent if Z is replaced by 7. 

At the one-loop order and in the limit of vanishing electron mass, the amplitude 
fully projects onto six invariants, in a similar way as in the THDM as follows: 

1=1 

where the invariants At are given by: 

1 ~t~ 75 1 — 75 

A = v{Pe+) APz)— U{pe^) , A2 = V{pe+) APz)— u{pe-) (2.1) 

1 -|- 1 — 
A = V{pe+) i>Z —U{p^-\pe-e{j)z)) , A = ^^(Pe+) ^ ^(Pe- ) (Pe" e(Pz) ) 

1 ~h 'T's 1 — '7^5 
A = ^^(Pe+) ^ u{Pe-){Pe+(^{Pz)) , A = ^^(Pe+) ^ u{Pe-){Pe+(^{Pz)) 

Here e is the polarization vector of the gauge boson Z. Summing over the Z gauge boson 
polarizations, the squared amplitude may be written as: 

4(M| - t)i?e[A1iA<; + M2MI] + (M| - t)'[|-M3|' + IM^l^] 
+4(M| - u)i?e[A^iA^; + M2M;] + (M| - ?i)'[|A<5|' + l-Mel'] 
-2(MlM| + M|s - tu)i?e[A<3-M5 + A^4A<6]} (2-2) 
The differential cross-section reads: 

^(e+e- - Z^O) = ^^^\/(^ - {Ma + MzY){s - (M^ - M,y) \M'\\ (2.3) 

In the on-shell scheme defined in it is found that there is no counter-term for the 
vertices Z-Z-A^, Z-7-A*^ and 7-7-74". Consequently the one-loop vertices Z-Z-A", Z-7-A'' 
and 7-7-^4° have to be separately UV finite. It can be seen from the expression of the 
invariant A (eqs. ^?T|) that the operators Ai^2 are dimension 3 while operators A,4,5,6 are 
dimension 5. One concludes that A^3,4,5,6 have to be UV finite while M.1^2 can potentially 
be UV divergent. We have checked numerically and analytically that the full set of 
invariants M.i are UV finite as required, and this feature will provide us with a good 
check of our calculation. 

3. Numerical results and discussion 

In this section we focus on the numerical analysis. We take the fine structure constant in 
the Thomson limit: a = 1/137.03598. 
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Figure 2: Total cross-section at 500 GeV centre of mass energy for tan /3 = 4 (left) and 
tan /3 = 40 (right) 



The MSSM Higgs sector is parametrized by the mass of CP-odd M4 and tan f3 while 
the top quark mass and the associated squark masses enter through radiative corrections 
2^ . In our study we will include the leading corrections only, where the light Higgs mass 



is given by: 



with 



1 r 



m 



AZ 



m\z — 4:m\mz cos^ (2/3) — 4e (m\ s| + m| 



+ m| + e 



3G 



F m 



Hog 



m| 

2 



(3.1) 



(3.2) 



The parameter m? = m^^m^^ denotes the average mass squared of the stop particles. A 
recent search for and A'^ excludes the region 0.5 < tan/3 < 2.4 [Q, assuming maximal 
mixing in the stop sector. In light of this we limit ourselves to the case where tan/3 > 3. 
Following the approach of [ES] we assume a detection threshold of 0.1 fb for e'^e~ J^Z . 
This would give 50 events before experimental cuts for the expected luminosities of 500 
fb~^. In the THDM this criterion would require tan/? < 0.3, even for a light pseudoscalar 
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The chargino/neutralino sector can be parametrized by the usual Mi, M2 and /x. We 
limit ourselves to the case of /i > since /i < seems to be disfavoured by 6 ^ S7 
constraint^ and from recent measurements of the anomalous magnetic moment of the 
muon, {g — 2)^ |3y]. We also use the SUGRA constraint Mi M2/2 and take the left 
selectron, right selectron and sneutrino to be degenerate with a common mass M-^ 

We present our results for two specific cases: (i) Light SUSY case (case L), where all 
sleptons, charginos and neutralinos are relatively light with mass < 200 GeV. (ii) Heavy 
SUSY case (case H) where all sleptons, charginos and neutralinos have mass > 450 GeV. 



^the chargino-stop contribution to 6 — > 57 add constructively with the SM one for /i < [E9[ . 
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• Case L: M2 = 220, = 180 and = 200 GeV 

• Case H: M2 = 500, /i = 600 and M; = 450 GeV 

We start by recalling that the THDM contribution to e+e^ A^Z in the small 
tan/3 regime is enhanced by the top quark contribution, leading to significant cross- 
sections (about a few fb). In the large tan/? regime the cross-section is suppressed and 
does not attain observable rates. In the MSSM we limit ourself to the case where tan (3 > 3, 
and consequently the THDM contribution is suppressed to the order of ~ 0.001 fb at 
\/s = 500 GeV. Our aim is to see if the SUSY contribution can enhance the cross-section. 

In Fig. 2 we show the dependence of the total cross-section on the CP-odd Higgs 
mass Ma for ^/s = 500 GeV in the case where tan/5 = 4 (Fig.S.a) and 40 (Fig.S.b). In 
both cases we show the MSSM contribution for case L and case H, as well as the THDM 
contribution. It can be seen that the light SUSY particle scenario (L) may enhance 
the total cross-section, which can reach ^ 0.003 fb for Ma < 300 GeV and tan/? = 4, 
with smaller cross-sections for tan/3 = 40. Note also that for Ma < 300 GeV the box 
contribution is an order of magnitude smaller than the vertex corrections in the MSSM 
case L, for both tan/5 = 4 and 40. The kink observed in this figure close to Ma ~ 2mt 
is a threshold effect due to the opening of the decay A° — »• ti. One can also have other 
threshold effects which arise from charginos and neutralinos. Note that in the case L and 
for tan/5 = 4 (resp 40) the mass of the light chargino is about 132.5 GeV (resp 146 GeV), 
and one can see a smooth threshold effect due to the opening of the channel — > xtxi 
for Ma ~ 264 GeV (resp Ma ~ 290 GeV). Due to the decoupling property of SUSY 
theory, the MSSM contributions in case H are very small and the total cross-section is 
close to the THDM one. 

We illustrate in Fig.3 the angular distribution in the MSSM with light SUSY particles 
for ^/s = 500 GeV, tan/? = 4 (left curve) and tan/5 = 40 (right curve). The shape of 
the graph is in agreement with that expected for a CP-odd scalar being produced by an 
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Figure 4: Total cross section a (fb) for e^e 
and MSSM case H (tan /3 = 4 right) 
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ZA, MSSM case L with tan /3 = 4 (left) 



effective ZZAP coupling [pT| . 

We stress here that in both cases, the shape of the angular distribution originates from 
the vertex contribution, since the box contribution is rather small at ^/s = 500 GeV. In 
the case of heavy SUSY particles (MSSM case H) the total cross-section is suppressed 
which would render difficult any analysis of the angular distribution. The numerical 
value of lies in the range 0.0019-0.0009 fb (resp 0.0008-0.0004 fb) with < ^ < tt 
for tan/? = 4 and Ma = 150 GeV (resp Ma = 250 GeV). We note in passing that the 
corresponding values for the THDM with tan (3 > 3 would also be very small. 

It is interesting to study the behaviour of the total cross-section versus i/i. Since a 
variety of SUSY particles are exchanged in the vertices and boxes, we expect that more 
threshold effects will appear. In Fig. 4, we present the total cross-section against y/s 
in the MSSM case L (left curve) and MSSM case H (right curve) for several values of 
Ma- It can be seen from the plot that there is a peak in both cases. In the MSSM case 
L, the peak appears at y/s ~ 400 GeV since in this case the masses of SUSY particles 
are approximately 200 GeV. In MSSM case H, the peak shows up for ^/s ~ 1 TeV. As 
one can see, the lighter Ma is, the more spectacular are the peaks. In the MSSM case 
L, for ^/s < 1.2 TeV the dominant contribution to the total cross-section is from vertex 
diagrams, while for energy ^/s > 1.2 TeV it is the box diagrams which dominate. Similarly 
in the MSSM case H, the vertex contibution dominates for ^/s < 2.2 TeV while the box 
diagram dominate if y/s > 2.2 TeV. 



4. Summary 

We have computed the cross-section for the production mechanism e~^e~ J^Z at high 
energy e^e~ colliders in the framework of the MSSM. Such a process proceeds via higher 
order diagrams and is strongly model dependent. 

The calculation was performed within the dimensional regularisation scheme and we 
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presented results for both the MSSM and THDM. In the former case hght SUSY particles 
may give important contributions to the cross-section, resulting in maximum values of 
order 0.003 fb. Therefore the SUSY enhancement is not sufficient to produce an observable 
signal at the planned luminosities of 500/6~^. In the MSSM with explicit CP violating 
phases, the pseudoscalar contains a CP even component and may be produced with an 
observable cross-section |32|. Therefore signals in the Higgsstrahlung channel for all h^, 
and could not be explained in the MSSM unless SUSY sources of CP violation 
are present. Observation of this process might enable one to distinguish between SUSY 
and non-SUSY Higgs sectors. We showed that threshold effects may enhance the cross- 
section. 
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